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ABSTRACT 

The binary system PSR B1259-63 consists of a 48 ms pulsar in a 3.4 year or- 
bit around a Be star and unpulsed TeV gamma-ray emission has been detected near 
previous periastron passages. A likely source of the gamma-rays is the inverse Comp- 
ton up-scattering of target photons from the Be star by the ultra-relativistic elec- 
tron/positron pulsar wind in the region of the plerion shock front. In this study the 
effect of the infrared emission from the Be star's circumstellar disc on inverse Compton 
gamma-ray production in PSR B1259-63 is investigated by considering an isotropic 
photon/electron approximation. The modified photon density distribution is deter- 
mined by using a curve of growth method fitted to previous optical and infrared 
observations. The inverse Compton scattering rate is calculated using the modified 
photon distribution and the exact scattering equation. It is shown that including the 
infrared emission increase the GeV gamma-ray flux by a factor > 2. 
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1 INTRODUCTION 

Be stars - defined as B-type stars showing or having showed 
emissi on lines inste ad of absorption lines in the Balmer spec- 
trum (|Collinslll987l ) - are known to produce an infrared (IR) 
fiux which is hi gher than pred icted by Kurucz stellar atmo- 
sphere models (|Kurucj|l979l ). It has long been speculated 
that Be sta rs possess an extended circumstellar envelope 
or disc (e.g. IStruveiri93ll ) and optical interferometry obser- 
vations confirm the presence of an extended circumstellar 
env elope which is symmetric al around the rotation axis (see 
e.g. iQuirrenbach et al.|[l994l ). These circumstellar discs act 
as reverse accretion discs as they are regions (confined to 
a disc structure) of increased stellar outflow. The free-free 
and free-bound emission which occurs in the disc is presum- 
ably responsible for the observed IR excess. These discs also 
show great variability and grow and shrink over periods of 
hundreds to t housands of days. For a r eview of classical Be 
stars, see e.g. IPorter fc Riviniud (|2003l ). 

An important Be X-ray Pulsar binary (Be-XPB) ex- 
ample is the TeV gamma-ray system PSR B1259-63 which 
consists of a 48 ms pulsar in orbit around a Be star, SS 2883, 
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and was initially detecte d as part of a radio survey of the 
southern Galactic plane (jjohnston et al.l Il992l . 1 19941 ) . The 
orbit is eccentric (e « 0.87) and unpulsed TeV gamma-ray 
emission has been detected by HES clo se to the perias- 
tron passage (Aharonian ct al. '20051. [2OO9I ). The circumstel- 
lar disc is mis-aligned to the orbital plane and the pulsar is 
believed to pass through it twice: before and after perias- 
tron. Previous observations have also detected unpulse d ra- 
dio emission near periastron (e.g lJohnston et al.lll999l ) and 
the system has been o bserved at X-ray frequencies across the 
whole orbit (see e.g. Cominskv. Roberts fc Johnston 1993 : 
iKaspi et al.ll 19951 : IChernvakova et al.ll2006l . l2009t ). These ob- 
servations show that the unpulsed emission is variable at all 
wavelengths with respect to the orbital period. 

The system is powered by the spin-down luminosity 
of the pulsar and the unpulsed radiation is believed to 
originate from a stand-off' shock front between the pulsar 
and stella r wind ; ext ensive modelling has been presented in 
[javani & AronsI |l997). In this model the X-rays are pro- 
duced via synchrotron radiation of ultra-relativistic elec- 



^ The High Energy Stereoscopic System, located near the Gams- 
berg Mountain in Namibia, consists of four Cherenkov imaging 



telescopes. 
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trons (Lorentz factor 7 ~ 10®) and the gamma-ray emis- 
sion through the inverse Compton (IC) scattering of target 
photons from the Be star. The ratio between the X-ray and 
gamma-ray flux is dependent on the radiative and adiabatic 
cooling times of the post-shocked wind and IC cooling was 
predicted to increase near periastro n. IC cooling of the pre- 
shocked wind was also suggested bv lBall fc Kiri3 1I2OOOI ). 

A slower wind model h as also been proposed 

jChernvakova fc lUariono^ Il999l . I2OO0I ) where the pulsar 
wind has a Lorentz factor of 7 ~ 10 — 100 and up-scatters 
the target photons from the Be star via IC processes to 
X-ray energies. The TeV gamma-rays are then explained 
via the interaction of the pulsar with the circumstellar 
disc and are the result of proton-prot on collisions, IC scat- 
tering and bremsstrahlu ng emission (|Kawachi et al.l |2004| : 
IChernvakova et aLlliooil ). 

If the TeV gamma-rays are the result of IC scattering 
of target photons from the Be star, this implies that elec- 
tron energies are of the order of TeV and the observed TeV 
gamma-rays are the result of scattering which occurred in 
the Klein-Nishina limit. While previous models have con- 
sidered IC scattering in the general case, there are currently 
only TeV gamma-ray observations of PSR B 1259-63 and the 
emphasis has been to model the TeV emission which will 
occur in the Klein-Nishina limit. Previous models have not 
considered how the scattering of the IR excess from the cir- 
cumstellar disc, occurring in the Thomson limit, will influ- 
ence the gamma-ray prod uction. 

It can be shown (e.g. iBlumenthal fc Gould|[l970l ) that 
IC scattering between a relativistic electron with a Lorentz 
factor 7 and a photon with a frequency v will occur in the 
classical Thomson limit if 



where h, nic and c represent Planck's constant, the elec- 
tron mass and the speed of light respectively. To explain 
the results of previous H ESS observations of PSR B1259-63, 
[Ah aronian et al ] liooi) proposed that electrons are acceler- 
ated in the shock between the pulsar (7 ~ 10® — 10^) and 
stellar wind, resulting in the IC scattering of stellar photons. 
This scattering will occur in the Thomson limit for target 
photons with frequency u <C 10^*(7/10®)^^ Hz, i.e. radio 
to IR. The IR excess provides an additional reservoir of soft 
target photons for IC scattering and an attractive framework 
to explain gamma-ray production through IC scattering in 
the Thomson limit, where the cross-section is significantly 
higher than the Klein-Nishina limit. A more detailed dis- 
cussion, especially focussing on the IC spectral properties, 
will be presented in section 3] 

Since the IR excess is tied to the growth and decay of 
the circumstellar disc, the size of the disc near periastron 
will influence gamma-ray production not only through the 
constraints on the plerion shock, but also through IC scat- 
tering in the Thomson limit. 

This study will discuss the effect of the IR excess on 
IC gamma-ray production in Be-XPBs using PSR B1259- 
63 as a trial case. Results of modelling the full IC scattering 
spectrum (Thomson and Klein-Nishina limit) of an isotropic 
distribution of photons and electrons will be presented here 
as an initial argument towards the importance of the con- 
tribution of the IR excess. The additional effects of the dif- 



ferent electron cooling processes in the plerion shock front 
and the changing scattering angle, due to orbital motion, 
will not be discussed at this time as the goal of this paper is 
to validate the importance of the IR excess on IC gamma- 
ray production. In order to include the IR excess the target 
photon distribution will be determined using the cur ve of 
growth method proposed by iLamers fc Waters! l|l984 ) (re- 
ferred to hereafter as the COG method) fitted to previous 
optical and IR observations. 

An outline of previous observations and studies of 
PSR B1259-63 will be presented (section[2]) before discussing 
the COG method (section [3]) and IC scattering (section Q. 
Section [5] will outline the modelling undertaken in this study 
and sections[6]and[7]will discuss the results and implications 
thereof. 



2 PSR B1259-63/SS 2883 

While this paper focuses only on one aspect of the the multi- 
wavelength emission from PSR B1259-63, a complete model 
needs to take into account all wavelengths to accurately ex- 
plain the radiative processes in the system. Below a brief 
summary is given of previous observations and modelling of 
the system which is relevant to the discussion presented in 
this paper. 

There have been few photometric observations of 
SS 2883, the optical companion to PSR B1259-63, but some 
measurements are available in bright star catalogues (Ta- 
ble[T]). For example. I Westerlund fc Garnie 3 (|l989l ') presented 
photometric and visual extinction measurements for SS 2883 
as part of their catalogue of OB stars near the South- 
ern Coalsack. Comparison of the catalogue's observations 
show little change in the optical magnit ude. From this and 
other observations [Johnston et al.l (| 19941 ) confined the spec- 
tral range of SS 2883 to a 09 - B2 type star and assumed 
a spectral type B2e, with M, ~ 10 M0 and i?* ~ 6 Hq. 
[Johnston et al.l ([l994[ ') also undertook spectroscopic obser- 
vations of the system and observed double peak H/?, H7, 
H(5 and He I 5876 emission lines. Analysis showed that the 
H/3 emission occurred at a radius of 8.5 R,, assuming a Ke- 
plerian circumstellar disc, which infers a circumstellar disc 
radius 7?disc > 8.5 i?,. 

There are also few observations in the mid-to- 
near IR, but observa tions are avai lable from the 2M ASS 
( Skrutskic ct al. 2006) and MSJtfl ([Price et al.[ [200 mis- 
sions. The IR data is shown in Table [H Due to the large 
uncertainty in the C-band measurement it was not consid- 
ered when modelling this system. 

In addition to the pulsed radio signal, an unpulsed 
radio signal has be en detected from th e system close to 
periastron (see e.g. [Johnston et al.[[2005l ). Observations of 
the system during consecutive periastron passages show an 
eclipse of the pulsar - and a decrease in the unpulsed syn- 
chrotron emission - roughly 20 days before and after pe- 
riastron passage as it passes through and behind the cir- 
cumstellar disc. For example, the pulsar eclipse lasted from 
16 days before until 18 days after periastron during the 

Data from 2MASS and tlie Midcourse Space Experiment 
(MSX) are available on-line through the NASA/IPAC Infrared 
Science Archive. 
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Table 2. IR data for PSR B1259-63 from 2MASS {J,H,Ks) and 
MSX (A,C). 



Band 


Wavelength 


Magnitude 


Flux 




fim 




Jy 


J 


1.235 


8.026 ±0.027 




H 


1.662 


7.699 ± 0.055 




Ks 


2.159 


7.248 ± 0.020 




A 


8.28 




0.2675 ±0.0136 


C 


12.13 




1.087 ±0.5859 



2000 passage (jConnors et al.ll2002h and from 15.8 days be- 
fore until 16.1 days after p eriastron during the 2004 pas- 
sage (| Johnston et al.ll2005l ). At 16 — 20 days from peri- 
astron the binary separation is « 40 — 50 R*, assuming 
Mpuisar = 1.4 Mq, which implies that the circumstellar 
disc must at least extend to these distances. While the syn- 
chrotron radio observations show many similarities there are 
marked differences in the flux levels during consecutive pe- 
riastron passages, attributed to variations in the circum- 
stellar disc and mat erial (see e.g. fig. 3 and discussion in 
I Johnston et~al]|2005l ). 

As a result of the alignment of the unpulsed radio, X- 
ray, and the 2004 HESS hght curves into two pealcs around 
the s ame two orbital phases (be fore and after periastron pas- 
sage) IChernvakova et"al] (|2006l ) suggested that the peaks in 
radio and X-ray, and the production of gamma-rays, was due 
to the pulsar passing through the disc. This implied a cir- 
cumstellar disc with a half-opening angle 18? 5, inclined 70° 
to the orbital plane. This suggested disc geometry is, how- 
ever, inconsistent with the radio eclipse (|Khangulvan et aU 
[2003). In addition, the detection of TeV gamma-rays 47 days 
before the 2007 periastron passage, mis-aligned to the or- 
bital phase position proposed above, suggests that IC scat- 
tering still plays a dom inant role in gamma-ray production 
l|Aharonian et al.ll2009l ). 

The exact position of the circumstellar disc and details 
behind the observed unpulsed radio, X-ray and gamma-ray 
emission (slow wind, 7 ^ 10^, or fast wind, 7 ~ 10*') remains 
unsolved and the observations of PSR B1259-63 by Fermi 
around the next periastron passage may help answer these 
questions. In this study we assume that the electrons have 
a Lorentz factor 7 ^ 10® and that the X-ray emission is 
the result of synchrotron radiation while the gamma-rays 
are produced through IC scattering. Under this assumption 
we present results of the full IC scattering of the Be star's 
photon spectrum taking into account the IR excess produced 
by the disc. 



3 THE LAMERS & WATERS CURVE OF 
GROWTH METHOD 

In this paper the COG method proposed by 
iLamers fc Water^ (|l984h has been chosen to model 
the IR excess from the circumstellar disc. This method 
allows for some prediction of the variability of the IR excess 
that will occur because o f the changing size of the d isc. The 
metho d, a s discussed in iLamers fc Waters! (|l984D . IWatersI 
l|l986l ) and lTelting et al.l ([19981 ). is summarised below. 

Under the COG method it is assumed that the circum- 



stellar disc has a half-opening angle extends to a radius 
7?disc, and has a power-law density profile that decreases 
with distance as 



P^r) = pa 



For a Be star-disc system viewed fa ce-on the flux ratio be- 
tween the disc and star is given by (|Telting et al.ll 19981 ^1 



-F';^,disc -Bi'(Tdisc) 



JR. 



[1- 



-2r,.{9)] 



2q dq, 



(2) 



where is the flux from the star and F^,disc is the flux 
from the disc, B,^{Tdisc) is the Planck blackbody function, 
Tdisc is the temperature of the disc, Iv,* is the flux from the 
appropriate Kurucz model for the star, and q is the impact 
parameter in units of the stellar radius. The optical depth, 
Tv{q), is given by (Waters 1986) 

r4g) = £;.,disc<z"'"+'C(n,e), 

where 



C{n,9) 



2n-2 7 

COS y dy, 



where 6, the half-opening angle, is in radians and i?i^,disc is 
defined as the optical depth parameter for the disc. This 
parameter is given by 

-Ei^.disc = XxXt,d 

where 

Xx = A2{(1 - e-'^''/*=^^'-)/(/i;./fcTdisc)} 
X {(7(1/, Tdi,c) + b(i^, Tdisc)} 

and 



4.923 X W^'^z^T-^/^fi-^^pl 



Rq 



where A is the wavelength, 1/ is the frequency, k is the Boltz- 
mann constant, g{u, Tdisc) + b{i^, Tdisc) is the sum of the free- 
free and free-bound gaunt factors|f| is the mean of the 
squared atomic charge, n the mean atomic weight (units of 
proton mass) and ^ the ratio of the number of electrons to 
the number of ions. 

The COG method can then be used to model a Be 
star's optical/IR flux by modifying a Kurucz model atmo- 
sphere at IR wavelengths with the flux ratio given by equa- 
tion (121), which is calculated by fitting the five parameters 
n, X,, /Jdisc, Tdisc and 9 to IR data. The Kurucz model 
atmosphere is chosen by a least-squares fit to the optical 
data, which is assumed to be less affected by the IR ex- 
cess. In practise, however, assumed values are taken for the 



^ In order to keep the notation consistent within this paper, Rdisc 
is in units of centimetres and the integration in equation ((2]| 
is between — (Rd iac/ R*)- This is different to equation (6) in 
iTelting et al. I lll998l) where the integration is between — -Rdisci 
but Rdi sc is in units of stella r radii. In addition, the optical depth, 
T{q), in lTelting et al. I l ll998t) (their equations (2) & (6)) IS equiv- 
alent to 2ti/ (q) as is given here, when the disc is viewed face on. 
* The free-free + free -bound gaunt factors published in tabular 
form in lWaters fc Lain crs ( 1984,) for wavelengths between 10~^ — 
6 cm, were used in this study. 
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Table 1. Optical magnitudes and colours for SS 2883 



Catalogue 


LSS 


Spectral Type 


V 


B-V 


U-B 


(U-B)o A, M,(U-B)o 


Wcstcrlund & Garnicr (1989) 


2883 


OB+ce,le,h 


10.01 


0.754 


-0.506 


-1.261 3.25 -4.6 


Klarc Sz Ncckcl C1977) 


2883 




10.04 


0.74 


-0.44 




Schild, Garrison & Hiltner ( 19831 


2883 




10.05 


0.72 


-0.47 




Drilling fl991~l 


2883 




10.07 


0.73 


-0.47 





disc temperature (taken to be a fraction of the stellar tem- 
perature) and half-openi ng an^le. The temperat ur e ratio, 
Tdisc = 0. 8 T,, used in iLamers fc Waters! (|l984l '). IWaterj 
l| 19861 ) andlTelting et all (Il998l ) has beenargued to be too 
high bv lMillar fc Marlboroughr(|l999h from their analysis of 
the energy gain/loss rate in four Be stars. For this reason 
a lower temperature ratio of Tdisc = 0.5 T, will be used 
in this study instead. From the pulsar eclipse a last con- 
straint can be placed on the disc radius and J?disc = 50 R, 
will be used in most instances. This corresponds to the 
separation between the pulsar and Be star at the start of 
the pulsar eclipse. In fact, from polarisation and rotation 
measurements made around the 1994 periastron passage, 
Ijohnston et al.l (|l996l ) suggested that the pulsar begins to 
pass into the disc ~ 100 days before and after periastron at 
a distance of ~ 150J?, . 

With this additional constraint on the disc radius there 
are only two free parameters {n & X,) which must be fitted 
to match the IR excess. 



4 INVERSE COMPTON SCATTERING 

A very brief discussion of IC scattering, presenting only 
the results relevant for t his study, will now be given . For 
detailed reviews see, e.g. B lumenthal fc Gouldl (|l970l ) and 
iRvbicki fc Lightman (2004). 

IC scattering involves the up-scatter of photons to 
higher energies by relativist ic electrons. It can be shown (e.g. 
iBlumenthal fc Gouldlll970l ) that the interaction between rel- 
ativistic electrons with velocity j3 = v/c and photons with 
energy hv results in the photons being up-scattered to en- 
ergies 

ft; A'y^hv 

in the Thomson limit {'yhv « rricC?), and 

2 

e~f fti jiricC 

in the extreme Klein-Nishina limit {-yhu » rricC^ and 
P — >■ 1). Although in the Thomson limit the characteris- 
tic energy of the scattered photon is increased by a factor 
^ 7^^ , it is still small compared to the energy of the relativis- 
tic electron. In the case of scattering in the ultra-relativistic 
limit the scattered photon energies are limited by the elec- 
tron energies. 

In a single scattering event (a single electron scattering 
a single photon) the energy of the electron before scatter- 
ing is given by '^m^c? and the energy of the photon, before 
and after scattering, by e and ei respectively. By expressing 
the energy of the scattered photon in units of the electron 
rest mass energy, i.e. E\ = ei/'ym^c? , the scattering rate per 
unit energy, for a single electron interacting with an isotropic 



photon density distribution n(e) (i.e. n umber density of pho- 
tons per unit energy), is expressed as (jBlumenthal fc Gouldl 
[19Z3) 



dN~/,e 2'KrQmeC n(e) de 



dtdEi 



2qlnq + {l + 2q){l-q) + l 



2 1-f r,g 



(3) 



where 



4e7 



El 



mec^' " r.(i-£;i)' 

and ro represents the classical electron radius. This is the 
exact expression for the scattering rate and is appropriate 
for all energies, provided that 7 S> 1, with Ei within the 
range 



In equation ((3|, n{e) de is the differential photon number 
density and in the case of a blackbody distribution of pho- 
tons is (e.g. IBlumenthal fc Gouldl[l970l ) 

1 ,2 

n(e) = 



To extend the calculations to multiple electron scatter- 
ings an electron distribution function needs to be considered. 
This distribution is normally expressed as a power law func- 
tion e.g. 



7Ve(7) = Kef 7mm < 7 < 7max, 

— 0, elsewhere, 



(4) 



where the distribution is confined to the energy region 7min 
to 7max. The total scattering rate per energy for an electron 
distribution is then found by intergrating over the initial 
electron and photon energies, i.e. 



dtdei 



dN-, 



'ynieC^ dtdEi 



df de. 



(5) 



5 MODELLING THE IC SCATTERING 

5.1 Modelling the Infrared Excess 

The optical and infrared dat a used to model PSR B 1 259-6 3 
in this paper is taken from IWesterlund fc Garnieil ^M), 
2MASS and MSX (Tables |T]fc|ll). The IR and optical data 
were de-reddened with the dipso s oftware package using 



the ex tinction measurement given in I Westerlund fc Garniei 
( 198S ) and following the assumptions of Ijohnston et al 



19941 ) 



The initial Kurucz atmosphere fits were done using 
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Table 3. Optical and IR fit to PSR B1259-63 



Optical fit 

Tcff 25000 K 

logg 3.5 

y^hift 19.64 

Curve of growth fit 

n 2.37 2.37 

logX. 7.87 7.32 

Rdisc 50 R* 50 R* 

Tdisc 12500 K 12500 K 

e 5° 18? 5 



the iKurucd (|l992t ) model atmospheres obtained from the 
Space Telescope Science Institut^f]. The tables only begin 
at a minimum frequency of ~ 3.3 x 10^'^ Hz and needed 
to be extended further into the IR regime. The last few 
data points lie within the Rayleigh- Jeans limit (hv -C kT) 
for a stellar temperature of ~ 20 000 K, and since in this 
limit log/i/,* oc \ogv, the data points were extended to 
V K. 5.2 X 10® Hz (or A ~ 6 cm) using a straight line fit. 

The stellar temperature, T*, and effective gravity, logg, 
were determined by finding the best least squares fit of the 
Kurucz model atmospheres to the optical data above 3 x 
10^* Hz (in order to limit the influence of the IR excess 
on the optical fit). The decrease in flux due to distance to 
the source from Earth (d) was also accounted for by scaling 
the intensity with th e parameter V^hift = log(7r_R^) — log(d^) 
l|Tehing et al.|[l998l ). A best flt temperature and effective 
gravity of T* = 25000 K and log g — 3.5, respectively, were 
found, consistent with an early B type star. 

The IR excess was then fltted to the IR data using 
the COG method. The disc temperature was chosen to be 
Tdisc = 0.5 T, = 12500 K and the disc radius _Rdisc = 50 R, 
as was inferred fro m the pulsar eclipse . In this paper the 
convention used by iTelting et al.l l|l998h was followed and 
the half-opening angle was held ai 9 = 5° . The remaining 
two parameters (n and X*) were simultaneously fitted with 
the Levenberg-Marquardt method using the implementation 
in the ORACLE software package. This produced a best fit of 
n = 2.37 and logX* = 7.87 and the resulting fit is plotted 
in Fig.[T] Fit s using the larger half-ope ning angle {9 ~ 18?5) 
suggested bv lChernvakova et al.l (|2006t l were also considered 
but produced no noticeable effect on the fitted optical/IR 
spectrum. The values are summarised in Table [3] 

As a futher check on these model simulations, COG 
fits were prod uce d using the da t a and parameters given in 
IWaterd (|l98d) and lTehing et al.l (|l998h for S Cen and X Per- 
sei respectively. These models were compared to the pub- 
lished results and were found to match extremely well. 

5.2 Modelling the IC gamma-ray spectrum 

The total number of IC scatterings is calculated from equa- 
tion ([5]) and the photon distribution 71(e) is derived from the 
predicted fiux as is given by the COG fit. This is different to 
previous models where either a blackbody or mono-energetic 
photon spectrum was assumed. The integration over e was 

^ http://www.stsci.edu/science/starburst/Kurucz.html 
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Figure 1. A Kurucz stellar atmosphere (solid fine) and the mod- 
ification create by the COG method (dashed line) fitted to data 
for PSR B1259-63/SS 2883. The open circles are the optical and 
IR data points. 

performed for the range e — 3.45 x 10~^^ —4.23 x 10"^^° ergs 
(~ 5.2 X 10^ - 6.4 X 10" Hz). At these limits, the ratio 
n(eii,„it)/max(n(e)) is 1.13 x 10"^ and 1.60 x 10"^® respec- 
tively. The integration over 7 was performed between the 
given 7min and 7max limits for each electron distribution 
case. 

The electron distribution, A'^e, in the IC calculation de- 
pends on where the electrons are assumed to originate, and 
what additional cooling processes are applicable. The exam- 
ples of pre-shocked, post-shocked adiabatic cooling and post- 
shocked radiative cooling are presented below. Since the aim 
of this study is to show the relative change in the flux of the 
gamma-ray spectrum, the models are not calibrated to fit 
the HESS data and two simplifying assumptions have been 
made in the analysis below; the value of = 1 was adopted 
for all cases and the photon density has not been scaled to 
the distance of the pulsar. Scaling the photon density sim- 
ply decreases the overall gamma-ray fiux by a constant value 
and does not change the relative shape. 

5.2.1 Pre-shocked electrons 

The possibility that IC scattering may occur in the pre- 
shocked region sur rounding the pulsar has been discussed by 
iBaU fc Kirk|||2000l) . In such a scenario the pulsar wind would 
have an approximate mono-energetic distribution around 
7 ^ 10^. Fig. [5] presents this scenario, which was calculated 
from equation © without integrating over 7. The scatter- 
ing rate shows a significant increase when the IR excess is 
included. 

5.2.2 Post-shock electrons - Adiabatic cooling 

If cooling is dominated by adiabatic cooling in the post- 
shock region the electron distribution is assumed to be 
a power law {Ne = 7~''). The HESS observations o f 
PSR B1259-63 suggest that p x 2.2 (|Aharonian et al.ll2005t ) 
and this was adopted as a first approximation. Two en- 
ergy ranges for 7 were considered, 7 — 10^ — 10^ and 
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^o'^ 10"'' 10"^ lo"'^ 10"' 10° 10"^ 10"* 10"^ io"^ io"' lo" 

Energy [TeV] Energy [TeV] 



Figure 2. The predicted total scattering rate, dNtot /dt dei, cal- 
culated for a mono-energetic electron 7 = 10^. The solid line is the 
expected gamma-ray flux without an IR excess and the dashed 
line shows the increased gamma-ray flux when the IR excess is 
included. 



Figure 3. The predicted total scattering rate, dNtot /dtdei, cal- 
culated using p = 2.2 for the range 7 = 10* — 10^. The electron 
distribution follow equation (|4]l. The solid line is the expected 
gamma-ray flux without an IR excess and the dashed line shows 
the increased gamma-ray flux when the IR, excess is included. 



7 = 10'' — 10^, encompassing a broad and narrow elec- 
tron distribution. The resulting scattering rates are shown 
in Fig. |3] and Fig. |4] respectively. For the broad electron 
distribution there is negligible influence on the gamma-ray 
flux when the IR excess is included. For the higher energy, 
narrow band, electron distribution, this is not the case and 
there is a significant increase in the scattering rate. 

A third example of a predominately adia- 
batic cooled electron distrib ution was taken from 
iKirk. Ball fc Skiaeraaseiil (|l999l ). Here the authors con- 
sidered p = 2.4 and 7 = 5.4 x lO'^ - 5.4 x lO'' and the 
sp ectrum is compar a ble to the HESS observations (see fig. 7 
m lAh aroman et al.l (Hooi)). The resulting scattering rate 
(assuming = 1) is shown in Fig. [5l There is a significant 
increase in the scattering rate due to the IR excess, though 
slightly less than is presented in Fig. [l] The peak in the 
gamma-ray spectrum also occurs at a lower energy. 

5.2.3 Post-shock electrons - Radiative Cooling 

When the electrons in the post-shocked region also cool via 
synchrotron radiation, the electron distribution can not be 
considered as a simple power law. While the post-shocked 
electrons initially have a power law distribution this is mod- 
ified as the electrons lose energy via synchrotron radiation. 
An example of dominant radiation cooling, with a magnetic 
field strength of B = 0.32 G, was considered, also taken 
ft-om iKirk et all (|l999l l. The initial electron spectrum was 
assumed to have an index of p = 1.4 and the energy range 
was 7 = 4.3 X 10^ — 4.3 x lO'^. This modified electron dis- 
tribution was used to model the scattering rate shown in 
Fig. [6l The details of the modification ar e beyond the scop e 
of this paper and the reader is referred to lKirk et al.l l|l999f ). 
As with the previous models a simplifying assumption of 
Ke = 1 has been made. The resulting plot is similar to 
Fig. |4] and Fig. [5l and shows a significant increase in the 
scattering rate. 

The electron distributions used are summarised in Ta- 
bled 



No disc 

50 R. disc 
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Figure 4. The predicted total scattering rate, dNtot /dtdei, cal- 
culated using p = 2.2 for the range 7 = 10'^ - 10''. The electron 
distribution follow equation (|4]l. The solid line is the expected 
gamma-ray flux without an IR excess and the dashed line shows 
the increased gamma-ray flux when the IR excess is included. 

6 DISCUSSION 

6.1 The Influence of the IR excess 

The results above show that for the case when higher energy 
electrons dominate the IC scattering, the gamma-ray spec- 
trum can be influenced by the infrared excess. Fig. [7] shows 
the increase in the gamma-ray flux due to the inclusion of 
the IR excess as the ratio Fu^star+disc/ Fu,star versus energy 
for the different electron distributions considered. 

The extent to which the gamma-ray spectrum is affected 
by the inclusion of the IR excess is dependent on which 
components of the Be star/disc spectrum are scattered in 
the Thomson limit. The peak in the stellar spectrum occurs 
at ~ (2 - 3) X lO^'^ Hz (Fig. [TJ, and IC scattering wiU 
occur in the Thomson limit for electrons with Lorentz factors 

7 < 10^ (iz/lO^'^Hz)"^ (see equation H}). For the case of 
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Table 4. Summary of electron distributions considered. 
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Mono-energetic 


pre-shocked 
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Broad distribution 


post-shock 


2.2 
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Narrow distribution 


post-shock 


2.2 
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10^ 
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Adiabatic (Kirk et al. 1999) 


post-shock 
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5.4 X 10'' 
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Radiative fKirk et al. 1999) 


post-shock 


1.4 
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Figure 5. The predicted total scattering rate, dNtot /dt dei, cal- 
culated using p = 2.4 for the range 7 = 5.4 X 10® — 5.4 X lO' (see 
iKirk et al.lll999l) . The electron distribution follow equation Q. 
The solid line is the expected gamma-ray flux without an IR 
excess and the dashed line shows the increased gamma-ray flux 
when the IR excess is included. 
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Figure 7. Fractional increase of the flux due to the inclusion 
for the IR excess for different electron distributions. The solid 
and dotted lines follow a power law distribution with p = 2.2, the 
dashed line is for a mono-energetic pulsar wind, and th e circle and 
crosse s correspond to the electron distributions from iKirk et al.l 
(1999) discussed in the text. A disc radius of 50 /?* is assumed 
for all cases. 



Energy [TeV] 

Figure 6. The predicted total scattering r ate, dNtot / dt dei, cal- 
culated using the electron spectrum from iKirk et al.l lll999l ) for 
dominate radiative cooling; B = 0.32 G and p = 1.4 for the range 
7 = 4.3 X 10^ —4.3 X 10'' . The solid line is the expected gamma-ray 
flux without an IR excess and the dashed line shows the increased 
gamma-ray flux when the IR excess is included. 



the broad electron energy distribution (7 — 10* — 10^) the 
gamma-ray spectrum is dominated by electrons with Lorentz 
factors 7 — 10* scattering photons from the peak of the 
stellar spectrum. This scattering will occur in the Thomson 
hmit and produce gamma-rays with energies 



3 X 1015 Hz 



eV. 



As a result, the major contribution at GeV energies is from 
the Thomson scattering of the whole stellar spectrum and 
the resulting gamma-ray spectrum is negligibly affected by 
the IR excess, as is shown in Fig. [3] and Fig. [T] 

For electrons with Lorentz factors 7 — 10^, the scat- 
tering will occur in the Thomson limit only at frequencies 
u <C 10^5 Hz. As a result, scattering of photons from the 
peak in the stellar spectrum will produce gamma-rays in the 
Klein-Nishina limit, while photons from the IR excess will 
scatter in the Thomson limit and increase the gamma-ray 
flux at GeV energies. 

In Fig. [7] it is shown that, with the exception of the 
broad (7 = 10® — 10^) electron distribution, the gamma-ray 
flux increases by a factor > 2 at energies less than a few GeV. 
The results show that the inclusion of the IR excess in the 
target photon distribution can influence the production of 
gamma-rays in PSR B1259-63, particularly at GeV energies. 
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Figure 8. Broad energy distribution of SS 2883 predicted using 
the COG method. The soUd line is the standard Kurucz model, 
while the dashed lines are the resulting IR excess for increasing 
disc sizes. The disc sizes are (from the bottom to the top) 10, 20, 
30, 40 and 50 -R». Only the disc radius is changed while the disc 
temperature and another parameters are the same as in Fig. [T] 



6.2 Variability of the Gamma-ray emission 

Since it is known that the circumstellar discs around Be stars 
are variable over periods of hundreds to thousands of days, 
the resulting change in the IR flux should also create vari- 
ability in the GeV gamma-ray emission from PSR B1259-63. 

A simple prediction can be made to show how the 
changing IR excess will modify the gamma-ray emission for 
PSR B1259-63. By changing the disc radius parameter in the 
COG method (while keeping the other parameters constant 
i.e. n, , Tdisc and 0) the magnitude of the IR excess can be 
simulated for different size circumstellar discs in PSR B1259- 
63. The resulting IR variability is shown in Fig. [S] where the 
lower solid line is a standard Kurucz model, and the dashed 
lines show the IR excess for increasing disc sizes. The con- 
sequential increase in the IC flux is shown in Fig. |9] as the 
ratio Fij^atar+disc/ Fv^star vcrsus cuergy The same disc sizes 
are used as in Fig.|S]and the electron distribution follows the 
power law used in Fig.|4] Fig. [9] shows that by including the 
IR excess in the target photon distribution, the gamma-ray 
flux can increase by a factor of > 2 and show variability in 
the GeV energy range. 

The modelling above assumes that the circumstellar 
disc only influences the photon distribution at IR wave- 
lengths. However, it is known that the circumstellar disc also 
increases the magnitude of Be stars at optical wavelengths. 
This will greatly increase the total number of target photons 
and may result in an even greater possible variability of the 
gamma-ray flux during consecutive periastron passages. An 
increase in the Be star's flux at optical wavelengths could 
also lead to variability in the case of the broad electron dis- 
tribution (7 = 10* - 10'') as an increase in the flux near the 
peak of the stellar spectrum will lead to an increase in the 
total gamma-ray flux. 
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Figure 9. Predicted increase in the IC flux for different disc radii, 
shown as the ratio of the modified flux to the initial flux. The disc 
sizes are 10, 20, 30, 40 and 50 -R» from the bottom to the top, 
and the photon number distribution is calculated as is shown in 
Fig. El 

6.3 Detectability of the gamma-ray modulation 

The increase in the gamma-ray flux due to the inclusion of 
the IR excess in an isotropic approximation will not be ob- 
servable with current gamma-ray telescopes. If the modelled 
fluxes are compared to the HESS spectrum the modified flux 
will fall below the Fermi detection threshold. A larger broad 
electron spectrum with more electrons around 7 ~ 10'' may 
produce a gamma-ray spectrum large enough to be be de- 
tected by Fermi, but the variability due to the IR excess 
reported here would not be observable. 

The analysis discussed here does however assume only 
isotropic scattering and the geometric size of the disc has not 
been considered. Anisotropic modelling of the system may 
show a more pronounced effect, especially close to periastron 
where the circumstellar disc not only provides an IR excess 
but scattering will occur over the larger solid angle of the 
disc. 



7 CONCLUSION 

PSR B1259-63 is a particularly important gamma-ray bi- 
nary system since it contains an independently confirmed 
48 ms pulsar compact object (as opposed to LS 5039 and 
LS I -|-61°303 where the nature of the compact object must 
be inferred) and the gamma-ray emission is most likely cre- 
ated through the IC scattering of target photons from the 
optical companion via the electrons/positrons emanating 
from the pulsar wind. Previous modelling of the IC gamma- 
ray emission has not taken into account the IR excess created 
by the circumstellar disc and how these additional target 

photons will infiuence the IC scattering. 

In this paper the iLamers fc Waters! (|l984l ') COG 
method was fitted to SS 2883 using optical, near-IR and mid- 
IR data available in literature and archives, under the as- 
sumption that the circumstellar disc has a radius of -Rdisc = 
50Rt , and that the Be star is a main sequence star with the 
parameters T, = 25000 K, log g = 3.5, M, ~ 10 M© and 
_R, ~ 6 R0 based on the least squares flt performed in this 
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paper and IJohnston et alj (| 19941 ). The fit was used to cal- 
culate the target photon number density which was in turn 
used to calculate the IC spectrum for PSR B1259-63. For a 
first approximation only isotropic scattering was considered 
in order to demonstrate the relevance of the IR excess. 

The modelling shows that the IR excess has a non- 
trivial influence on the IC spectrum, particularly at GeV 
energy gamma-rays where the flux increases by a factor > 2 
for reasonable electron distributions. It has also been shown 
(since the circumstellar disc is variable) how changes in the 
IR excess will manifest as variability at gamma-ray ener- 
gies. Consequently, observations of consecutive periastron 
passages may show differences in the GeV gamma-ray light 
curve as a result of a possibly changing IR excess. 

The results presented above have implications for other 
Be-XPBs and gamma-ray binary systems. Previous mod- 
elling of the IC spectrum only considered a blackbody or 
mono-energy distribution for the target photons, however, 
the results above show that the IR excess from the cir- 
cumstellar disc will have a non-trivial influence on the IC 
gamma-ray spectrum that must be accounted for when mod- 
elling PSR B1259-63. 
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